SUSY Flat Directions - to get a VEV or not? 
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Abstract. We investigate the potential of SUSY flat directions (FDs). Large FD vacuum expecta- 
tion values (VEVs) can delay thermalisation and solve the gravitino problem - if FDs decay per- 
turbatively. This depends on how many and which directions get the VEVs. Recently the decay of 
the FDs have been studied with the VEVs as input. Here we look at how the VEVs come about - 
statistically and analytically. 
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Supersymmetry, MSSM and flatness 

Supersymmetry, in which all ferminic particles get a bosonic partner and vice versa, 
is a very well motivated extension of the Standard Model (SM). It gives unification 
of all the gauge couplings at the same energy scale and it automatically removes all 
the high energy divergencies from Quantum Field Theory due to the equal magni- 
tude, opposite sign contributions of fermionic and bosonic loops in the high energy 
limit. No superpartners have been observed. This emberrassment can be hidden by 
introducing R-parity, where all the known particles (including Higgs bosons) get +1, 
and all superpartners -1. Thus superpartners can only be created and destructed in 
pairs. This immediately makes the lightest superpartner (LSP) the favorite Dark Mat- 
ter candidate - indeed, giving the LSP typical weak scale mass and coupling gives 
an energy density in the right ballpark. The scalar potential of the MSSM (the mini- 
mal supersymmetric extension of SM) consists of F-terms (^^ \ dW /d^\'^) and D-terms 

(Lfl,?fl/2| ^^r"^ p) where T", W are the gauge generators and the Superpotential re- 
spectively - and are the scalar fields. The renormali sable part of the superpotential is 
[1] Wrenorm = yj UiQjHu - y'J DiQjHd - yj EiLjHd + IiHuH/ where the Yukawa cou- 
plings are the same as in the Standard ModeP. 

Flatness means that the potential is zero for nonzero field values - which can be seen to 
happen if, and only if, all D-terms and F-terms vanish individually. This can only happen 
for exact (unbroken) SUSY and with no nonrenormalisable terms in the superpotential. 



Superfields and their scalar part will be represented by the same symbol. Q,L,E,U,D,Hu,Hj are 
lefthanded quarks, lefthanded leptons, righthanded charged leptons, up-type righthanded quarks, down- 
type righthanded quarks, positive-hypercharged Higgs, negative-hypercharged Higgs. 
^ We will choose the basis where the SUSY-breaking mass terms - not the Yukawas - are diagonal. 



Evolution of flat directions 



A catalogue of FDs where given in [2]. FD evolution was studied in [3]. Giving VEVs^ 

Q\'=Ll = Dl=(j>/V3 (1) 

The superfields can be multiplied and the product can be parameterised by a canonical 
field that experiences Hubble friction 

X - QI'lIdI X = c(l>'^ (m = 3), ^ + 3H^+V'{(j>) = (2) 

m is used to keep example as general as possible. Adding nonrenormalisable terms to 
the superpotential (M is a breaking scale: Planck/GUT/other) 

/L 

^ — ^renorm ~l~ j^/n— 3 ^ 

where all possible gauge invariant and R-parity conserving terms will be allowed (and 
expected to be of order 1). All FDs can be lifted by such terms - either by itself 
■;^X - ;^^<^>" (n = m - if positive R-parity), itself squared -j^Z^ = ;^^<^>" 
{n — 2m - if negative R-parity) or a combination of fields in the flat direction with exactly 
one field not in the direction -j^^i^c^"^ (with respect to which the derivative then can 
be taken). The potential is: 

V(0)=m^,|0p-ctf^|0p+(^ ^ '/^^^J + (4) 

The masses are the "real" masses, whereas the terms of the same order in (j) are Hubble 
induced terms. The F-term (last term) has order 2n — 2 because it is the derivative squared 
of n'th order (W„ contributes to V2n-2)- The A-term (third term) is a coupling only 
between scalars. During inflation the induced terms dominate. If c > (50% chance?? - 
however, with minimal Kaehler potential it does NOT work), the minimum is displaced 

to (^^^^ — ^ " ^ (j8 order 1 constant) which could easily be of order lO^^GeV, while 

after inflation the real masses dominate and the VEV will oscillate around a minimum 
of zero. The phase dependence of the potential creates Baryon number if 5 — L is not 
conserved by the FD. 

Sometimes flatness is described by monomials. The relation is this: L1L2E2 = (Vg * 
jU — e * V^) * T*^ is a monomial. This gives flatness to either term (Ve, jU, t'^) or (e, V^, t'^) 
which can be ^(e^*^\e'*^,e'*^^) which keeps Da = for all generators. Superterm 
W4 oc L\L2Et,Ni gives \Fn^ P \L\L2E'i\^ which is positive and thus lifts flatness, whereas 
the A-term is A * e'*^^LiL2E2, + h.c. which is negative for one of the mentioned field 
combinations and thus chooses the minimum. 



^ Flatness: hypercharge sum to zero, one up and one down weak charge balance, and a color and same 
anticolor balance - lower index: generation, upper index: color, weak charge (in that order for Q). 



Cosmological consequences 



[4] pointed out that FDs induce masses to inflaton decay products of order g| {^) \ >Hj 
(g: gauge coupling, Hi: Hubble parameter during inflation). This prevents preheating"^ 
since MSSM scalars are no longer massless. This lowers the cosmological reheating 
temperature (from 10^ to 10^ — lO^GeV) and this avoids the gravitino problem (why we 
don't see any). Even, the FD could be the inflaton itself[5]! However, [6] noticed that 
FD is only important if it lives long enough i.e. decay perturbatively - which it claimed 
was unlikely. Rather there would be immediate decay of the FD through quick particle 
production. 



Framework and particle production 

In [7] we developed the framework for analysing particle production from FDs. The 
crucial thing is to work in the unitary gauge, where no unphysical Goldstone bosons 
appear. 

Writing the excitations of fields in a vector: S = {t,\...E,i...E,nY , we can write the 
Lagrangian as 

^|a^S|2-^S^^22_s^[/S + ... (5) 

Make an orthogonal transformation S' = AS (A is orthogonal) with 

A^A = U (6) 
and the mixed kinetic term (U-term) disappears: 

^^lld^Z'l^-^Z'^^'^Z' (7) 

where = AJ^'^A^ = ABJ^jB^A^ = C^^C^ , and C = AB. 
Non-perturbative particle production is investigated like this [8]: 
Change to conformal fields Xi = ci'^'p where a denotes the scale factor with equation 

of motion 

Xi + ^jj{t)Xj = (8) 
where dots represent derivatives with respect to conformal time t, and 

ajj = a^^'l + k^5ij (9) 

where k labels the comoving momentum. Using an orthogonal time-dependent matrix 
C{t), we can diagonalise via {t)Q?'{t)C{t) = a)^(?),where (D is diagonal. 



^ The instantanious transfer of energy from the inflaton to massless scalars through a parametric reso- 
nance. 



As the vacuum changes, a new set of creation/annihilation operators are required. We 
use Bogolyubov transformation with Bogolyubov coefficients a and j8 (matrices when 
more than one field). 

Initially a = I and j8 = while the system evolves as (matrix multiplication implied): 

a = -i(oa + -la -jp 
2(0 

$ = —a + i(oP-Ja-ip, (10) 
2(0 

with the matrices / and J given by 

2 \ ^ ^/a J 

J = 1 ( x/roC^C^-^C^C^y (11) 

2 \ yUO J 

Occupation number of the ith bosonic eigenstate reads (no summation) 

mit) = (j8*i3^),-. (12) 

This showes (still [8]) that not just rapidly changing eigenvalues of the Hamiltonian, but 

also rapidly changing eigenstates can create particles from the vacuum. 

Since initially a = I and /3 = 0, eq.lO shows that a non- vanishing matrix i is a 
necessary condition to obtain /3 7^ and hence ^ 0. In our framework we have^ 

C^C = B^A^AB = -B^UB (13) 

which shows that we just need to find U and ^ and diagonalise the latter - with no need 
to make explicit transformation to the primed system. 



A single FD 

One flat direction often mentioned in the literature is LLE. 
The potential is 

y = i (^Djj + Y,Dl^ withDH = Y^qim^ andDa = y^^P'(I> (14) 

where qi is the hypercharge, and g 1,^2 are the hypercharge- and SU2 gauge couplings. 
We give these VEV's: 

{li) = (pe'''^ (Vm)=0 (15) 



The last equation only holds if B is constant in time. 



The Lagrangian reads 

^=i\ \D,^i\' kv - 1 W^y (16) 

i=l i 

with F, W are hypercharge- and weak field strenth tensors and where for field 0j: Df = 
{d^ — iqiA^)5ij — La=i ^Pfj^a 0j is the covariant derivative. is the Pauli-matrix. 
We get mixed kinetic terms^ 

-(p%((7i+C72-2a3)-(p2A3(di-C72) (17) 

- these (diagonal) gauge- VEV derivative mixtures are unphysical Nambu-Goldstone 
bosons (Goldstones). Making a f/(l) gauge transformation on the multiplets (4>,) 



4>, = e^i^'^, with A = 



(18) 



and by making a 5?/ (2) gauge transformation 



«D,. ^ «d;. = with 7 = ^ (19) 



we get rid of the Goldstones and the VEVs are 



in) = (pe'"" (20) 



with (7 = ((7i + <72 + <73)/3. We write excitations 



jU = ((p + ^3)/'"'^^ y^^ = {^7 + me''' (21) 

Among the kinetic terms we find 

^D-cp (Ai (4 + 4) +A2(6 - 15)) • (22) 



" The gauge fields have a suppressed Lorentz index: zero - only this component matters, if the phases are 
assumed to change rapidly in time, not in space. 



- again we have Goldstones - this time from the off-diagonal gauge generators. They are 
removed by redefinitions 

v. = (,p + &)/"^' e=(^iM."' 

v2 

= + v,= «5^.'- (23) 

and we are in the unitary gauge! We calculate U, diagonalise the latter and find 7 = 

- no particle production. 

We found in [7] in a toy model that particle production is propertional to the deriva- 
tive of phase differences between the participating VEV fields. Here we gauged both 
differences away and found no particle production. 

UDD (< >= (^e'^i, < >= < >= (j)e''^^) works exactly as LLE. Both 
phase differences are gauged away - no particle production. In {QQQ)^\L2L2,E^ . 
(Squarks with identical 5?/ (2) -charge chosen for simplicity) the VEV-fields are (can 
be chosen to be) 

(<p + ^io)/^'^^^ fi = ((p + ^g)/^''"^''''"*^''''"^^ (24) 
and the no- VEV fields (notice complicated normalisation needed) 



^ 4: transform as a 4 under SU (2) i.e. uncontracted SU (2) indices. 



,,c2 _ |ii+i|i2 „(<Ti „cl _ "all i^\2 Jay 

■K, — ( fog+i'i^ao _ lai+ffeA J(cti+ct2+ct3) 
~ I, 2^/5 V30 j ^ 



;,c3 _ ( fe7+'fe8 I fog- ''^30 I fol-'fe \ „('(Ti /TCX 

' 2,/I 2V5 ^ ^/30 I ^ '^^^^ 



27+'S28 I <,29-lh30 i S31-'S32 | ^IGi 
2V3 2\/5 \/30 





_ ^fol+ifo2 


d'^ = 


I ^/2 




( ^\5~i^l6 1 < 
I V2 ' 



i~'fe2 I ^n-i^ii feg+'feo foi+'fo2 \ i(a\—2as) 

V6 2V3 2V5 V30 y 



V3(-^/fe+^/^+4^) , 

An example of nonzero entry in 7 is ^4 17 = 73 ig = ^ ^ . The structure 

4\/l0(fe2+6g2,p2)5 



is general in the relavant (p » k limit: n'x ^JgijG- (n: particle density, k: momentum 
of produced particle) - a huge number proportional to a VEV phase difference. 



Several FDs 

UDD and LLE can coexist. Combined there are 6 VEV fields, and one can only gauge 
4 phase differences away (4 diagonal generators). However, it is just one phase for each 
direction - and they don't interact^. There is no particle production. But LLE and QLD 



U is block diagonal and / = 0. 



can also coexist with an L field in common. 



(26) 



where A is the relation between the magnitude of VEVs. This gives particle production. 



There are problems with this picture. Monomials (or directions) are not independent. 
There are only 17 mass terms - or 20 if righthanded neutrinos(A^s) are included - yet 
there are 712 (715 including A'^s) independent monomials [10](counted after my talk). 



independent. Also, if{QQQ)4LiL2Li,Ei has a VEV, so has {QQQ)4LiL2L2Ei. Also, 
when is QQQLLLE broken? From earlier arguments one could imagine that without 
A^s it would be broken by itself squared - dimension 14. However, the space of Q, L, E is 
27(18+6+3) dimensional. It breaks the Standard Model completely, so D-terms remove 
12 complex degrees of freedom (c.d.o.f.)'^. So the D-flat space is 15 dimensional. W4 (4th 
order superpotential) includes QQQL and QULE - so Fq,Fl,Fu ,Fe give 36 complex 
constraints and thus W4 lifts the flat direction. This means it is lifted by the 6th order 
in the potential- eventhough its A-term is of much higher order. Including A'^s will give 
A-terms like QQQLLLE N but the direction will still be lifted by W4 (including LLEN). 
So the relation between flat directions and monomials has really broken down - see [10]. 



The potential must be investigated for the following reasons. It is very well to state that 
particle production is proportional to VEV phase differences. But do these differences 

have dynamical equations of motion to drive them? Also, the effective mass term must 
be negative for any direction to get a large VEV"'. There are 712 monomials - but 
also combinations thereof are gauge invariant {LLE, UDD but also LLE * UDD can 
get couplings). The formally flattest direction (the one lifted by highest order in W) is a 
combination of Q, U, E (from monomials UUUEE, QQQQU, QUQUE). It is only lifted 
by W9 (Vi6)[2]. Include A^s, and it is lifted by We (Vio). (Just add N to the mentioned 



Problems with this picture 



To illustrate, m^^ 

2 

^{QQQULiL2L2Ei 




Investigation of the potential 



One real non-flat direction and one real gauge choice for each. 
' This can be avoided by large A-term [11]. 



monomials [10].) Also, while [5] claims that VEVs are in general hierarchical or flat 
directions independent, [6] claims that there will be several large ones. So we aspire to 
write down the general potential to 10th order (a rough count: 2.3 million couplings). 



Normalisation and statistical treatment 

We want to count the number of couplings correctly. The 2^^ made of 3 Qs look 
like (GGG)"^ = QfQfQl"eabc£pY " it is 8 (not 27) dimensional [2]. For instance 
QQQm + QQQni + QQQin = QQQm - QQQin = while QQQn^+QQQni-iQQQm 
is a free parameter. If one let all 3 Qij^k with two I's and one 2 be standard gaussian 
N{0, 1) they will have N{0, 1 / \/3) projected in the relevant direction. Adding the three 
will get us back to A'^(0, 1) - but we do want to know that there are one, not three, 
parameters - and it is only the same distribution when gaussianity is assumed. [2] has 28 
types of monomials (without family indices). We have found that these can be combined 
to ~700 gauge invariant combinations of less than or exactly 10 fields - in ~400 unique 
field combinations. (Field combination HuHjLLE have combinations HuHj + LLE and 
HuL + LH^E. The dimension even of a monomial can be nontrivial. {QQQ)2 combined 
with QU has neither 8*9=72 (product of dimensions) nor 12*3=36 (12 ways to assign 
at least 2 different generations to 4 Qs and 3 generations ofU- but rather 54 dimentions 
(stated in [2], written down in [10]). For normalisation we will choose the antisymmetric 
tensors of SU (2), 5(7 (3) -contractions to have norm 1 - while other linear combinations 
of field vectors - ie. combination of family indices - will be treated as if they were 
basis vectors (ie. if a^b are products (including 5[/(2),Sf/(3)-contractions) of fields, 
{a — b)/ Vl will be used as basis vector if a + b = 0). Finally, we've chosen l/n\ for 
each superfield appering n times in a product (must be superfield, since not welldefined 
for fields i.e. (HuHa)^ = {H+R-f + IH+R-H^H^ + [H^H^f cannot be normaUsed 
by 1 and 1 /2 simultaniously). 



Work to do, in progress 

There is a statistical approach: Choose random couplings. Find minimum of potential 
with Monte Carlo methods. Try enough combinations to get a feeling of what VEVstruc- 
ture is typical. There is an analytical approach: Impose symmetries. Assume common 
couplings (mi/2,ino,A and so on...). Investigate the role of the (formally) flattest direc- 
tion. I work currently on both approaches. 



" Doublet under ^f/ (2). 



Conclusion 



SUSY Flat Directions can have crucial influence on (p)reheating and offer a very nice 
solution to the gravitino problem, baryogenesis and even offer a "known" particle as a 
candidate for being the inflation. Preheating is a serious threat to this. The jury is still 
out, and the potential must be investigated thoroughly. 
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